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Abstract. Precise efforts in theoretical astrophysics are needed to fully 
understand the mechanisms that govern the structure, stability, dynam- 
ics, formation, and evolution of differentially rotating stars. Direct com- 
putation of the physical attributes of a star can be facilitated by the use 
of highly compact azimuthal and separation angle Fourier formulations 
of the Green's functions for the linear partial differential equations of 
mathematical physics. 



1. Introduction 

Solutions to the three dimensional inhomogeneous linear partial differential 
equations of mathematical physics are expressible in terms their appropriate 
Green's functions (Duffy 2001). Single, double, and triple eigenfunction inte- 
gration and summation expressions for Green's functions are obtainable in the 
coordinate systems that allow for separation of variables (Morse & Feshbach 
1953). Separation of variables maps the linear three dimensional homogeneous 
partial differential equations into three decoupled ordinary differential equations 
(Miller 1977). Important linear 3D partial differential equations of mathematical 
physics include the Laplace equation, the Helmholtz equation, the biharmonic 
equation, the wave equation, the diffusion equation, and the Schroedinger equa- 
tion. The properties of R-separability for these partial differential equations 
are reflected in the Green's function expansions. In the rotationally invariant 
coordinate systems that R-separably solve these equations, new discrete Fourier 
representations exist for the Green's functions for the Laplace equation (Cohl et 
al. 2000). Discrete Fourier expansions, given in terms of the azimuthal and sep- 
aration angles, must exist for the rest of the linear partial differential equations 
of mathematical physics. The coefficients of these discrete Fourier representa- 
tions will be given in terms of identifiable transcendental functions obtained by 
reversing and collapsing Green's function expansions for these partial differential 
equations. The fundamental mathematical tools required in order to complete 
this investigation are commonly available in the mathematics and physics liter- 
ature. 
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The Green's function for Laplace's equation is 

gL(x ' x ' ) = R^r (1) 

where x and x' are given by the Cartesian coordinates (x,y,z) and (x',y',z r ). 
Traditionally, the Green's function is expanded in terms of spherical harmonics 

g L = P^os 7 ), (2) 

in spherical (r,(f>,6) coordinates, where r< (r>) is the smaller (larger) of the 
distances r and r', cos 7 = cos cos 9' + sin0sin0'cos((/> — 4>'), and Pi is the 
Legendre polynomial of the first kind. 1 

By utilizing Heine's toroidal identity (Cohl & Tohline 1999, Cohl et al. 
2000, Cohl et al. 2001): 
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where |£| > 1, and, < ip < 2tt, one can derive two highly compact Fourier 
series representations for the Green's function of Laplace's equation 
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in spherical and cylindrical (R, (ft, z) coordinates respectively. Similar discrete 
Fourier transform expressions must exist for the linear partial differential equa- 
tions of mathematical physics. Algorithmic implementation of compact Fourier 
representations for the Green's functions of mathematical physics significantly 
improves the accuracy of the solution for these partial differential equations. 



2. The Biharmonic, Triharmonic, and Higher Harmonic Equations 

One may also express potential problems in terms of an inhomogeneous bihar- 
monic equation whose source is proportional to the Laplacian of the density. The 
solution is expressible in terms of an integral of it's Green's function, |x — x'|, 
convolved with source. The resulting inhomogeneous biharmonic equation 

V 4 $(x) = 4vr V 2 p(x), (6) 



1 "...The set of six coordinates {x, x'} may be viewed either as defining two points relative to the 
origin or as the coordinates of a three-body system once the center of mass has been separated. 
In the body frame, three out of the six coordinates are dynamical variables, the potential energy 
depending on them." (Cohl et al. 2001) 
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where V 2 is the Laplacian, represents the solution to the potential problem over 
the infinite three-dimensional domain. The Green's function for the biharmonic 
equation (Vautherin 1972) 

S 2L (x,x') = ±|x-x'|, (7) 

is a much better behaved integration kernel than for Laplace's equation because 
the singularity has been moved to infinity through 

$(x) = ~\J V' 2 p(x') |x-x'|. (8) 

One can generalize Heine's reciprocal square root identity using (Abramowitz & 
Stegun, 1965; Gradshteyn & Ryzhik 1994-to within a sign) 



r • • • r Qm-dowr = n tty - i) n/2 ^-i(o, o) 

J( J( 2 T(m + n + |J 2 

and by integrating both sides of eq. (2) with respect to £. Performing the integra- 
tions over £ we obtain the following general result for any integer n - generalizing 
the last expression at the bottom of p. 182 in Magnus, Oberhettinger, & Soni 
(1966) (see also eqs. (26) & (30) in Cohl et al. 2000) 

(C-cos^ = #^4 E T{m - n+ \\ Q n l(C )e^, (10) 
which reduces for n = 1 to 

7^^=^— e H-i) ^-i(o^. (ii) 
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Using this expression, we can rewrite the Green's function for the biharmonic 
equation in two different ways 



2<5 2L = |x-x'| = ^ T i Y (n 2 --) 1 Q 1 J r ^±ll\e^ (12) 

and 



2fe = * t (» ! l)"9l-i(^^)«^, (13) 
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This new formulation is amenable in analytic and computational physics appli- 
cations on high-performance-computing architectures, since Laplacians can be 
readily computed on computational fluid dynamics meshes. Higher order har- 
monic Green's function expansions can be generated using this method. Bound- 
ary values can now be computed effectively along an arbitrarily chosen z-axis 
(m = 0) and at values chosen to lie within the outer-most extent of a chosen vol- 
umetric region (m > 0). Physical solution of the interior inhomogeneous prob- 
lem can be obtained by proper boundary value treatment as described above. 
Accelerations can be obtained from the potential by performing a precise nu- 
merical gradient. In the case of cylindrical and spherical coordinates, the three 
dimensional Poisson solve can be further facilitated through the use of a dis- 
crete azimuthal Fourier transform. This decouples the three dimensional Pois- 
son problem into a set of decoupled 2D problems that can be solved with second 
order accurate finite differencing using either direct or iterative methods. 



3. Applications 

Using highly compact representations of infinite domain Green's functions, solu- 
tions to the linear partial differential equations of mathematical physics can be 
more easily obtained. In this paper, and in recent papers, we treat the Green's 
function for the Laplace equation. Here we have extended this result to the 
three dimensional biharmonic and higher harmonic equations. Further variants 
are possible for the 3D Helmholtz, wave, and diffusion equations. In the future, 
we intend to investigate and more precisely describe these new Green's function 
expansions. Many areas of theoretical physics will benefit greatly from precise 
numerical implementations of these compact expansions. 

Azimuthal Fourier identities for potentials such as Coulomb (|x — x'| _1 ) or 
Yukawa ( e ~ fc l x ^ x l| x — x'| _1 ) lead to new classical and quantum energy theo- 
rems. These allow for rapid and precise evaluation of the Coulomb or Yukawa 
direct and exchange interactions (Cohl et al. 2001). Classical Yukawa eigenfunc- 
tions are obtained through transcendental function identification of the Lamb- 
Sommerfeld integral (see Magnus et al. 1966). In quantum physics this is ac- 
complished through the use of the azimuthal selection rule for the self-energies, 
namely for the direct and exchange Hamiltonian elements only the m = and 
m = mi — nT-2 terms survive respectively. The application of the selection rule 
allows for exact evaluation of the Hamiltonian matrix elements for two-electron 
interactions in atomic physics, molecular physics, condensed matter physics, 
physical chemistry, and biology. Two-electron interactions are critically impor- 
tant in obtaining opacities and correct equations of state in astrophysically dense 
atomic and molecular fluid media. Magnetohydro dynamic problems can now be 
easily handled with Heine expansions for the Green's functions of potential the- 
ory. By expressing the equations of fluid dynamics of a compressible media in 
terms of a velocity potential and a vector Poisson equation, one may compute 
precise velocity boundary values in vortex and shock flow regions (Lamb 1932) 
by solving the appropriate Poisson problem. Radiation transport, and classical 
and quantum scattering, will be greatly facilitated through compact Poisson for- 
mulations of the Green's function for the 3D diffusion and Helmholtz equations. 
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Precise Coulomb and Yukawa energies in the nuclear Hamiltonian allow for a 
higher degree of precision in obtaining nuclear structure. 



4. Conclusion 



The ultimate resolution of these analytical investigations will be efficient algo- 
rithmic implementations of these new schemes. We propose these methods to 
the three dimensional star community in hope that you may continue to enjoy 
significantly improved economical and precise boundary values for studies of 
analytical and numerical three dimensional stellar astrophysics. 
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A Expressions for Unit Order Toroidal Functions 

The unit order negative one half degree Legendre function of the second kind 
can be expressed as 

Q -^-7W^) E (^Ti)- (15) 

where E is the complete elliptic integral of the second kind. The unit order, 
positive one half degree Legendre function can be expressed as 

where K is the complete elliptic integral of the first kind. We can express this 
same function in terms of the complete elliptic integral D 



Ql(0 = - 




(17) 
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as well. Higher degree, unit order toroidal functions of the second kind can be 
easily derived using the following recurrence relation, 

«^i<« = 2^1 «™-i«> "Sri «L|(0- 
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